Astrophysics in the Far Ultraviolet 

ASP Conference Series, Vol. **VOLUME**, YEAR 

Sonneborn, Moos & Andersson 



Semi-empiric Radiative Transfer Modeling of FUSE 
Stellar Spectra 



A. Lobel 1 , E. H. Avrett 1 , and J. P. Aufdenberg 2 

1 Harvard- Smithsonian Center for Astrophysics, 60 Garden Street, 
Cambridge 02138 MA, USA 

2 National Optical Astronomy Observatory, P.O. Box 26732, Tucson, 
85726 AZ, USA 

Abstract. We present an overview of radiative transfer modeling efforts to 
interpret spectra of a variety of stellar objects observed with FUSE. Detailed 
radiative transfer modeling of high ion emission line profiles of C in and O vi 
observed in the far-UV spectrum, provides a powerful means to probe the ther- 
mal and dynamic properties of high-temperature plasmas in the atmospheres 
of stars. We model asymmetric emission lines of Cm A977 (and Mgn h & k) 
observed in spectra of luminous cool stars such as a Aqr, to infer the wind- and 
microturbulence velocity structures of the upper chromosphere. Semi-empiric 
radiative transfer models that include transition region temperature conditions, 
are further developed based on detailed fits to O vi resonance emission lines in 
the supergiant a Aqr, the classical Cepheid variable (3 Dor, and to self-absorbed 
O vi emission lines in the cataclysmic variable SW UMa. 

We observe that the C III resonance line profile of a Aqr assumes a remark- 
able asymmetric shape, reminiscent of P Cygni type profiles observed in hot 
luminous supergiants. The model calculations indicate outflow velocities above 
~140 kms" 1 at kinetic temperatures of 65 kK and higher. Based on detailed 
model fits to the narrow red-shifted and self-absorbed O vi emission lines of SW 
UMa we compute that the gas- and electron density exceed the density conditions 
of the upper solar transition region by about three orders of magnitude. We pro- 
pose that the large gas density of p~1.4 10~ n grcm~ 3 favors a region of warm 
dense plasma of 100 kK <T gas < 300 kK that collapses onto the white dwarf with 
a mass accretion rate of 1—2 10 15 gs" 1 above or between the accretion disk. We 
discuss how detailed semi-empiric fits to emission lines observed with the high 
spectral resolution of FUSE can provide reliable constraints on the mass loss 
or mass accretion rates in these objects. The complete text of this conference 
presentation is available at cfa-www. harvard. edu/^alobel/conference. html 



1. FUSE Observations and Detailed Radiative Transfer Modeling 

We present a comparative study of transition region (TR) wind dynamics in a 
Aqr (G2 lb), SW UMa (CV), Dor (F-G la-lab) and the Sun, to develop semi- 
empiric radiative transfer models of the thermal structure of the stellar chromo- 
sphere and TR, and to determine the velocity and electron density structures. 
FUSE spectra of the hybrid supergiant a Aqr have been observed for Science 
Team Program P218 (Dupree et al. 2005). The classical 5-Cepheid variable j3 
Dor was observed in Aug. and Oct. 2003 for GI-D107 (PI A. Lobel), while the 
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Figure 1. Radiative transfer best fit (solid line) to C in A977 in a Aqr. 



dwarf nova SW UMa was observed in Nov. 2001 for GI-B074. The left-hand 
panel of Fig. 1 shows non-LTE radiative transfer fits to C in A977 through a 
semi-empiric model of the chromosphere and lower TR of a Aqr. The outward 
decrease of p (panels right) produces a scattering core in the computed emission 
profile (Lobel & Dupree 2000; 2001), which assumes an asymmetric shape due to 
opacity in a fast accelerating wind. Best multi-level atom model fits are obtained 
for V w i n d>140 kms -1 in the lower TR (T~65 kK), which strongly scatters the 
blue emission line wing. The model signals a supersonic optically thick warm 
wind in the outer atmosphere of this hybrid supergiant (Lobel & Dupree 2002). 

The left-hand panels of Fig. 2 show O vi lines in the cataclysmic variable 
SW UMa. The lines are remarkably far red-shifted with respect to heliocen- 
tric rest. We suggest that the double-peaked narrow profiles result from self- 
absorbed emission line formation in an inf ailing region near the white dwarf that 




Figure 2. Best fit (long-dash dotted line) to O vi A1032 & A1037 in SW UMa. 
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Figure 3. Best fit to solar O vi (upper left panel) emission compared to (3 Dor. 



is sufficiently optically thick at T gas between 100 kK and 300 kK (panels right). 
The depth of the O vi A1032 central absorption line core can only be computed 
with electron densities N e at least 1100 times larger compared to the solar TR 
values. The downflow velocity (Pwmd) in the model is computed from the conser- 
vation of total mass in spherical geometry with M = 4 7r p(r) V w i n d(r) (r+i?wD) 2 , 
using a mass accretion rate M of 1 10 15 gs -1 , and i?wD=0.01 R©. This M- value 
provides the best fit to the wavelength position of the central self-absorption core 
where V^i n d=+70 to +65 kms -1 in the O vi lines formation region. 

The upper left-hand panel of Fig. 3 compares the O vi lines we compute 
with an atmospheric model of the average quiet Sun with SOHO-SUMER ob- 
servations (solid dots). The lines form over a rather small region of ~10 model 
layers with 100 kK < T gas < 300 kK, or ~200 km into the solar TR above the up- 
per chromosphere (panels right). A best fit to the FWHM and equivalent width 
of the lines is obtained with projected microturbulence velocities (T4urb) increas- 
ing from at least ~30 kins -1 to ~39 kms -1 over this thermal range (Warren 
et al. 1997). A best fit to the core position of the O vi A1032 line is computed 
with a mean downflow velocity of +5 kms -1 . Preliminary atmosphere models 
for (3 Dor (lower left panel) require V tuv ^ exceeding 40 kms -1 in the O vi lines 
formation region to match the large line widths we also observe in a Aqr. 
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